3.7.16 \(\int x^2 (d+e x^2)^3 (a+b \arcsin (c x)) \, dx\) [616]

3.7.16.1 Optimal result
3.7.16.2 Mathematica [A] (verified)
3.7.16.3 Rubi [A] (verified)
3.7.16.4 Maple [A] (verified)
3.7.16.5 Fricas [A] (verification not implemented)
3.7.16.6 Sympy [A] (verification not implemented)
3.7.16.7 Maxima [A] (verification not implemented)
3.7.16.8 Giac [B] (verification not implemented)
3.7.16.9 Mupad [F(-1)]

3.7.16.1 Optimal result

Integrand size = 21, antiderivative size = 287 \[ \int x^2 \left (d+e x^2\right )^3 (a+b \arcsin (c x)) \, dx=\frac {b \left (105 c^6 d^3+189 c^4 d^2 e+135 c^2 d e^2+35 e^3\right ) \sqrt {1-c^2 x^2}}{315 c^9}-\frac {b \left (105 c^6 d^3+378 c^4 d^2 e+405 c^2 d e^2+140 e^3\right ) \left (1-c^2 x^2\right )^{3/2}}{945 c^9}+\frac {b e \left (63 c^4 d^2+135 c^2 d e+70 e^2\right ) \left (1-c^2 x^2\right )^{5/2}}{525 c^9}-\frac {b e^2 \left (27 c^2 d+28 e\right ) \left (1-c^2 x^2\right )^{7/2}}{441 c^9}+\frac {b e^3 \left (1-c^2 x^2\right )^{9/2}}{81 c^9}+\frac {1}{3} d^3 x^3 (a+b \arcsin (c x))+\frac {3}{5} d^2 e x^5 (a+b \arcsin (c x))+\frac {3}{7} d e^2 x^7 (a+b \arcsin (c x))+\frac {1}{9} e^3 x^9 (a+b \arcsin (c x)) \]

output
-1/945*b*(105*c^6*d^3+378*c^4*d^2*e+405*c^2*d*e^2+140*e^3)*(-c^2*x^2+1)^(3 
/2)/c^9+1/525*b*e*(63*c^4*d^2+135*c^2*d*e+70*e^2)*(-c^2*x^2+1)^(5/2)/c^9-1 
/441*b*e^2*(27*c^2*d+28*e)*(-c^2*x^2+1)^(7/2)/c^9+1/81*b*e^3*(-c^2*x^2+1)^ 
(9/2)/c^9+1/3*d^3*x^3*(a+b*arcsin(c*x))+3/5*d^2*e*x^5*(a+b*arcsin(c*x))+3/ 
7*d*e^2*x^7*(a+b*arcsin(c*x))+1/9*e^3*x^9*(a+b*arcsin(c*x))+1/315*b*(105*c 
^6*d^3+189*c^4*d^2*e+135*c^2*d*e^2+35*e^3)*(-c^2*x^2+1)^(1/2)/c^9
 
3.7.16.2 Mathematica [A] (verified)

Time = 0.16 (sec) , antiderivative size = 231, normalized size of antiderivative = 0.80 \[ \int x^2 \left (d+e x^2\right )^3 (a+b \arcsin (c x)) \, dx=\frac {315 a x^3 \left (105 d^3+189 d^2 e x^2+135 d e^2 x^4+35 e^3 x^6\right )+\frac {b \sqrt {1-c^2 x^2} \left (4480 e^3+80 c^2 e^2 \left (243 d+28 e x^2\right )+24 c^4 e \left (1323 d^2+405 d e x^2+70 e^2 x^4\right )+2 c^6 \left (11025 d^3+7938 d^2 e x^2+3645 d e^2 x^4+700 e^3 x^6\right )+c^8 \left (11025 d^3 x^2+11907 d^2 e x^4+6075 d e^2 x^6+1225 e^3 x^8\right )\right )}{c^9}+315 b x^3 \left (105 d^3+189 d^2 e x^2+135 d e^2 x^4+35 e^3 x^6\right ) \arcsin (c x)}{99225} \]

input
Integrate[x^2*(d + e*x^2)^3*(a + b*ArcSin[c*x]),x]
 
output
(315*a*x^3*(105*d^3 + 189*d^2*e*x^2 + 135*d*e^2*x^4 + 35*e^3*x^6) + (b*Sqr 
t[1 - c^2*x^2]*(4480*e^3 + 80*c^2*e^2*(243*d + 28*e*x^2) + 24*c^4*e*(1323* 
d^2 + 405*d*e*x^2 + 70*e^2*x^4) + 2*c^6*(11025*d^3 + 7938*d^2*e*x^2 + 3645 
*d*e^2*x^4 + 700*e^3*x^6) + c^8*(11025*d^3*x^2 + 11907*d^2*e*x^4 + 6075*d* 
e^2*x^6 + 1225*e^3*x^8)))/c^9 + 315*b*x^3*(105*d^3 + 189*d^2*e*x^2 + 135*d 
*e^2*x^4 + 35*e^3*x^6)*ArcSin[c*x])/99225
 
3.7.16.3 Rubi [A] (verified)

Time = 0.66 (sec) , antiderivative size = 287, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.238, Rules used = {5230, 27, 2331, 2123, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int x^2 \left (d+e x^2\right )^3 (a+b \arcsin (c x)) \, dx\)

\(\Big \downarrow \) 5230

\(\displaystyle -b c \int \frac {x^3 \left (35 e^3 x^6+135 d e^2 x^4+189 d^2 e x^2+105 d^3\right )}{315 \sqrt {1-c^2 x^2}}dx+\frac {1}{3} d^3 x^3 (a+b \arcsin (c x))+\frac {3}{5} d^2 e x^5 (a+b \arcsin (c x))+\frac {3}{7} d e^2 x^7 (a+b \arcsin (c x))+\frac {1}{9} e^3 x^9 (a+b \arcsin (c x))\)

\(\Big \downarrow \) 27

\(\displaystyle -\frac {1}{315} b c \int \frac {x^3 \left (35 e^3 x^6+135 d e^2 x^4+189 d^2 e x^2+105 d^3\right )}{\sqrt {1-c^2 x^2}}dx+\frac {1}{3} d^3 x^3 (a+b \arcsin (c x))+\frac {3}{5} d^2 e x^5 (a+b \arcsin (c x))+\frac {3}{7} d e^2 x^7 (a+b \arcsin (c x))+\frac {1}{9} e^3 x^9 (a+b \arcsin (c x))\)

\(\Big \downarrow \) 2331

\(\displaystyle -\frac {1}{630} b c \int \frac {x^2 \left (35 e^3 x^6+135 d e^2 x^4+189 d^2 e x^2+105 d^3\right )}{\sqrt {1-c^2 x^2}}dx^2+\frac {1}{3} d^3 x^3 (a+b \arcsin (c x))+\frac {3}{5} d^2 e x^5 (a+b \arcsin (c x))+\frac {3}{7} d e^2 x^7 (a+b \arcsin (c x))+\frac {1}{9} e^3 x^9 (a+b \arcsin (c x))\)

\(\Big \downarrow \) 2123

\(\displaystyle -\frac {1}{630} b c \int \left (\frac {35 e^3 \left (1-c^2 x^2\right )^{7/2}}{c^8}-\frac {5 e^2 \left (27 d c^2+28 e\right ) \left (1-c^2 x^2\right )^{5/2}}{c^8}+\frac {3 e \left (63 d^2 c^4+135 d e c^2+70 e^2\right ) \left (1-c^2 x^2\right )^{3/2}}{c^8}+\frac {\left (-105 d^3 c^6-378 d^2 e c^4-405 d e^2 c^2-140 e^3\right ) \sqrt {1-c^2 x^2}}{c^8}+\frac {105 d^3 c^6+189 d^2 e c^4+135 d e^2 c^2+35 e^3}{c^8 \sqrt {1-c^2 x^2}}\right )dx^2+\frac {1}{3} d^3 x^3 (a+b \arcsin (c x))+\frac {3}{5} d^2 e x^5 (a+b \arcsin (c x))+\frac {3}{7} d e^2 x^7 (a+b \arcsin (c x))+\frac {1}{9} e^3 x^9 (a+b \arcsin (c x))\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {1}{3} d^3 x^3 (a+b \arcsin (c x))+\frac {3}{5} d^2 e x^5 (a+b \arcsin (c x))+\frac {3}{7} d e^2 x^7 (a+b \arcsin (c x))+\frac {1}{9} e^3 x^9 (a+b \arcsin (c x))-\frac {1}{630} b c \left (\frac {10 e^2 \left (1-c^2 x^2\right )^{7/2} \left (27 c^2 d+28 e\right )}{7 c^{10}}-\frac {70 e^3 \left (1-c^2 x^2\right )^{9/2}}{9 c^{10}}-\frac {6 e \left (1-c^2 x^2\right )^{5/2} \left (63 c^4 d^2+135 c^2 d e+70 e^2\right )}{5 c^{10}}+\frac {2 \left (1-c^2 x^2\right )^{3/2} \left (105 c^6 d^3+378 c^4 d^2 e+405 c^2 d e^2+140 e^3\right )}{3 c^{10}}-\frac {2 \sqrt {1-c^2 x^2} \left (105 c^6 d^3+189 c^4 d^2 e+135 c^2 d e^2+35 e^3\right )}{c^{10}}\right )\)

input
Int[x^2*(d + e*x^2)^3*(a + b*ArcSin[c*x]),x]
 
output
-1/630*(b*c*((-2*(105*c^6*d^3 + 189*c^4*d^2*e + 135*c^2*d*e^2 + 35*e^3)*Sq 
rt[1 - c^2*x^2])/c^10 + (2*(105*c^6*d^3 + 378*c^4*d^2*e + 405*c^2*d*e^2 + 
140*e^3)*(1 - c^2*x^2)^(3/2))/(3*c^10) - (6*e*(63*c^4*d^2 + 135*c^2*d*e + 
70*e^2)*(1 - c^2*x^2)^(5/2))/(5*c^10) + (10*e^2*(27*c^2*d + 28*e)*(1 - c^2 
*x^2)^(7/2))/(7*c^10) - (70*e^3*(1 - c^2*x^2)^(9/2))/(9*c^10))) + (d^3*x^3 
*(a + b*ArcSin[c*x]))/3 + (3*d^2*e*x^5*(a + b*ArcSin[c*x]))/5 + (3*d*e^2*x 
^7*(a + b*ArcSin[c*x]))/7 + (e^3*x^9*(a + b*ArcSin[c*x]))/9
 

3.7.16.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 2123
Int[(Px_)*((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] 
:> Int[ExpandIntegrand[Px*(a + b*x)^m*(c + d*x)^n, x], x] /; FreeQ[{a, b, c 
, d, m, n}, x] && PolyQ[Px, x] && (IntegersQ[m, n] || IGtQ[m, -2])
 

rule 2331
Int[(Pq_)*(x_)^(m_.)*((a_) + (b_.)*(x_)^2)^(p_.), x_Symbol] :> Simp[1/2   S 
ubst[Int[x^((m - 1)/2)*SubstFor[x^2, Pq, x]*(a + b*x)^p, x], x, x^2], x] /; 
 FreeQ[{a, b, p}, x] && PolyQ[Pq, x^2] && IntegerQ[(m - 1)/2]
 

rule 5230
Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))*((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_ 
)^2)^(p_.), x_Symbol] :> With[{u = IntHide[(f*x)^m*(d + e*x^2)^p, x]}, Simp 
[(a + b*ArcSin[c*x])   u, x] - Simp[b*c   Int[SimplifyIntegrand[u/Sqrt[1 - 
c^2*x^2], x], x], x]] /; FreeQ[{a, b, c, d, e, f, m}, x] && NeQ[c^2*d + e, 
0] && IntegerQ[p] && (GtQ[p, 0] || (IGtQ[(m - 1)/2, 0] && LeQ[m + p, 0]))
 
3.7.16.4 Maple [A] (verified)

Time = 0.27 (sec) , antiderivative size = 404, normalized size of antiderivative = 1.41

method result size
parts \(a \left (\frac {1}{9} e^{3} x^{9}+\frac {3}{7} d \,e^{2} x^{7}+\frac {3}{5} d^{2} e \,x^{5}+\frac {1}{3} d^{3} x^{3}\right )+\frac {b \left (\frac {c^{3} \arcsin \left (c x \right ) e^{3} x^{9}}{9}+\frac {3 c^{3} \arcsin \left (c x \right ) d \,e^{2} x^{7}}{7}+\frac {3 c^{3} \arcsin \left (c x \right ) d^{2} e \,x^{5}}{5}+\frac {\arcsin \left (c x \right ) c^{3} x^{3} d^{3}}{3}-\frac {35 e^{3} \left (-\frac {c^{8} x^{8} \sqrt {-c^{2} x^{2}+1}}{9}-\frac {8 c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{63}-\frac {16 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{105}-\frac {64 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{315}-\frac {128 \sqrt {-c^{2} x^{2}+1}}{315}\right )+105 d^{3} c^{6} \left (-\frac {c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{3}-\frac {2 \sqrt {-c^{2} x^{2}+1}}{3}\right )+135 d \,c^{2} e^{2} \left (-\frac {c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{7}-\frac {6 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {8 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {16 \sqrt {-c^{2} x^{2}+1}}{35}\right )+189 d^{2} c^{4} e \left (-\frac {c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{5}-\frac {4 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{15}-\frac {8 \sqrt {-c^{2} x^{2}+1}}{15}\right )}{315 c^{6}}\right )}{c^{3}}\) \(404\)
derivativedivides \(\frac {\frac {a \left (\frac {1}{3} d^{3} c^{9} x^{3}+\frac {3}{5} d^{2} c^{9} e \,x^{5}+\frac {3}{7} d \,c^{9} e^{2} x^{7}+\frac {1}{9} e^{3} c^{9} x^{9}\right )}{c^{6}}+\frac {b \left (\frac {\arcsin \left (c x \right ) d^{3} c^{9} x^{3}}{3}+\frac {3 \arcsin \left (c x \right ) d^{2} c^{9} e \,x^{5}}{5}+\frac {3 \arcsin \left (c x \right ) d \,c^{9} e^{2} x^{7}}{7}+\frac {\arcsin \left (c x \right ) e^{3} c^{9} x^{9}}{9}-\frac {e^{3} \left (-\frac {c^{8} x^{8} \sqrt {-c^{2} x^{2}+1}}{9}-\frac {8 c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{63}-\frac {16 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{105}-\frac {64 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{315}-\frac {128 \sqrt {-c^{2} x^{2}+1}}{315}\right )}{9}-\frac {d^{3} c^{6} \left (-\frac {c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{3}-\frac {2 \sqrt {-c^{2} x^{2}+1}}{3}\right )}{3}-\frac {3 d^{2} c^{4} e \left (-\frac {c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{5}-\frac {4 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{15}-\frac {8 \sqrt {-c^{2} x^{2}+1}}{15}\right )}{5}-\frac {3 d \,c^{2} e^{2} \left (-\frac {c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{7}-\frac {6 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {8 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {16 \sqrt {-c^{2} x^{2}+1}}{35}\right )}{7}\right )}{c^{6}}}{c^{3}}\) \(417\)
default \(\frac {\frac {a \left (\frac {1}{3} d^{3} c^{9} x^{3}+\frac {3}{5} d^{2} c^{9} e \,x^{5}+\frac {3}{7} d \,c^{9} e^{2} x^{7}+\frac {1}{9} e^{3} c^{9} x^{9}\right )}{c^{6}}+\frac {b \left (\frac {\arcsin \left (c x \right ) d^{3} c^{9} x^{3}}{3}+\frac {3 \arcsin \left (c x \right ) d^{2} c^{9} e \,x^{5}}{5}+\frac {3 \arcsin \left (c x \right ) d \,c^{9} e^{2} x^{7}}{7}+\frac {\arcsin \left (c x \right ) e^{3} c^{9} x^{9}}{9}-\frac {e^{3} \left (-\frac {c^{8} x^{8} \sqrt {-c^{2} x^{2}+1}}{9}-\frac {8 c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{63}-\frac {16 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{105}-\frac {64 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{315}-\frac {128 \sqrt {-c^{2} x^{2}+1}}{315}\right )}{9}-\frac {d^{3} c^{6} \left (-\frac {c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{3}-\frac {2 \sqrt {-c^{2} x^{2}+1}}{3}\right )}{3}-\frac {3 d^{2} c^{4} e \left (-\frac {c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{5}-\frac {4 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{15}-\frac {8 \sqrt {-c^{2} x^{2}+1}}{15}\right )}{5}-\frac {3 d \,c^{2} e^{2} \left (-\frac {c^{6} x^{6} \sqrt {-c^{2} x^{2}+1}}{7}-\frac {6 c^{4} x^{4} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {8 c^{2} x^{2} \sqrt {-c^{2} x^{2}+1}}{35}-\frac {16 \sqrt {-c^{2} x^{2}+1}}{35}\right )}{7}\right )}{c^{6}}}{c^{3}}\) \(417\)

input
int(x^2*(e*x^2+d)^3*(a+b*arcsin(c*x)),x,method=_RETURNVERBOSE)
 
output
a*(1/9*e^3*x^9+3/7*d*e^2*x^7+3/5*d^2*e*x^5+1/3*d^3*x^3)+b/c^3*(1/9*c^3*arc 
sin(c*x)*e^3*x^9+3/7*c^3*arcsin(c*x)*d*e^2*x^7+3/5*c^3*arcsin(c*x)*d^2*e*x 
^5+1/3*arcsin(c*x)*c^3*x^3*d^3-1/315/c^6*(35*e^3*(-1/9*c^8*x^8*(-c^2*x^2+1 
)^(1/2)-8/63*c^6*x^6*(-c^2*x^2+1)^(1/2)-16/105*c^4*x^4*(-c^2*x^2+1)^(1/2)- 
64/315*c^2*x^2*(-c^2*x^2+1)^(1/2)-128/315*(-c^2*x^2+1)^(1/2))+105*d^3*c^6* 
(-1/3*c^2*x^2*(-c^2*x^2+1)^(1/2)-2/3*(-c^2*x^2+1)^(1/2))+135*d*c^2*e^2*(-1 
/7*c^6*x^6*(-c^2*x^2+1)^(1/2)-6/35*c^4*x^4*(-c^2*x^2+1)^(1/2)-8/35*c^2*x^2 
*(-c^2*x^2+1)^(1/2)-16/35*(-c^2*x^2+1)^(1/2))+189*d^2*c^4*e*(-1/5*c^4*x^4* 
(-c^2*x^2+1)^(1/2)-4/15*c^2*x^2*(-c^2*x^2+1)^(1/2)-8/15*(-c^2*x^2+1)^(1/2) 
)))
 
3.7.16.5 Fricas [A] (verification not implemented)

Time = 0.26 (sec) , antiderivative size = 277, normalized size of antiderivative = 0.97 \[ \int x^2 \left (d+e x^2\right )^3 (a+b \arcsin (c x)) \, dx=\frac {11025 \, a c^{9} e^{3} x^{9} + 42525 \, a c^{9} d e^{2} x^{7} + 59535 \, a c^{9} d^{2} e x^{5} + 33075 \, a c^{9} d^{3} x^{3} + 315 \, {\left (35 \, b c^{9} e^{3} x^{9} + 135 \, b c^{9} d e^{2} x^{7} + 189 \, b c^{9} d^{2} e x^{5} + 105 \, b c^{9} d^{3} x^{3}\right )} \arcsin \left (c x\right ) + {\left (1225 \, b c^{8} e^{3} x^{8} + 22050 \, b c^{6} d^{3} + 31752 \, b c^{4} d^{2} e + 25 \, {\left (243 \, b c^{8} d e^{2} + 56 \, b c^{6} e^{3}\right )} x^{6} + 19440 \, b c^{2} d e^{2} + 3 \, {\left (3969 \, b c^{8} d^{2} e + 2430 \, b c^{6} d e^{2} + 560 \, b c^{4} e^{3}\right )} x^{4} + 4480 \, b e^{3} + {\left (11025 \, b c^{8} d^{3} + 15876 \, b c^{6} d^{2} e + 9720 \, b c^{4} d e^{2} + 2240 \, b c^{2} e^{3}\right )} x^{2}\right )} \sqrt {-c^{2} x^{2} + 1}}{99225 \, c^{9}} \]

input
integrate(x^2*(e*x^2+d)^3*(a+b*arcsin(c*x)),x, algorithm="fricas")
 
output
1/99225*(11025*a*c^9*e^3*x^9 + 42525*a*c^9*d*e^2*x^7 + 59535*a*c^9*d^2*e*x 
^5 + 33075*a*c^9*d^3*x^3 + 315*(35*b*c^9*e^3*x^9 + 135*b*c^9*d*e^2*x^7 + 1 
89*b*c^9*d^2*e*x^5 + 105*b*c^9*d^3*x^3)*arcsin(c*x) + (1225*b*c^8*e^3*x^8 
+ 22050*b*c^6*d^3 + 31752*b*c^4*d^2*e + 25*(243*b*c^8*d*e^2 + 56*b*c^6*e^3 
)*x^6 + 19440*b*c^2*d*e^2 + 3*(3969*b*c^8*d^2*e + 2430*b*c^6*d*e^2 + 560*b 
*c^4*e^3)*x^4 + 4480*b*e^3 + (11025*b*c^8*d^3 + 15876*b*c^6*d^2*e + 9720*b 
*c^4*d*e^2 + 2240*b*c^2*e^3)*x^2)*sqrt(-c^2*x^2 + 1))/c^9
 
3.7.16.6 Sympy [A] (verification not implemented)

Time = 1.23 (sec) , antiderivative size = 525, normalized size of antiderivative = 1.83 \[ \int x^2 \left (d+e x^2\right )^3 (a+b \arcsin (c x)) \, dx=\begin {cases} \frac {a d^{3} x^{3}}{3} + \frac {3 a d^{2} e x^{5}}{5} + \frac {3 a d e^{2} x^{7}}{7} + \frac {a e^{3} x^{9}}{9} + \frac {b d^{3} x^{3} \operatorname {asin}{\left (c x \right )}}{3} + \frac {3 b d^{2} e x^{5} \operatorname {asin}{\left (c x \right )}}{5} + \frac {3 b d e^{2} x^{7} \operatorname {asin}{\left (c x \right )}}{7} + \frac {b e^{3} x^{9} \operatorname {asin}{\left (c x \right )}}{9} + \frac {b d^{3} x^{2} \sqrt {- c^{2} x^{2} + 1}}{9 c} + \frac {3 b d^{2} e x^{4} \sqrt {- c^{2} x^{2} + 1}}{25 c} + \frac {3 b d e^{2} x^{6} \sqrt {- c^{2} x^{2} + 1}}{49 c} + \frac {b e^{3} x^{8} \sqrt {- c^{2} x^{2} + 1}}{81 c} + \frac {2 b d^{3} \sqrt {- c^{2} x^{2} + 1}}{9 c^{3}} + \frac {4 b d^{2} e x^{2} \sqrt {- c^{2} x^{2} + 1}}{25 c^{3}} + \frac {18 b d e^{2} x^{4} \sqrt {- c^{2} x^{2} + 1}}{245 c^{3}} + \frac {8 b e^{3} x^{6} \sqrt {- c^{2} x^{2} + 1}}{567 c^{3}} + \frac {8 b d^{2} e \sqrt {- c^{2} x^{2} + 1}}{25 c^{5}} + \frac {24 b d e^{2} x^{2} \sqrt {- c^{2} x^{2} + 1}}{245 c^{5}} + \frac {16 b e^{3} x^{4} \sqrt {- c^{2} x^{2} + 1}}{945 c^{5}} + \frac {48 b d e^{2} \sqrt {- c^{2} x^{2} + 1}}{245 c^{7}} + \frac {64 b e^{3} x^{2} \sqrt {- c^{2} x^{2} + 1}}{2835 c^{7}} + \frac {128 b e^{3} \sqrt {- c^{2} x^{2} + 1}}{2835 c^{9}} & \text {for}\: c \neq 0 \\a \left (\frac {d^{3} x^{3}}{3} + \frac {3 d^{2} e x^{5}}{5} + \frac {3 d e^{2} x^{7}}{7} + \frac {e^{3} x^{9}}{9}\right ) & \text {otherwise} \end {cases} \]

input
integrate(x**2*(e*x**2+d)**3*(a+b*asin(c*x)),x)
 
output
Piecewise((a*d**3*x**3/3 + 3*a*d**2*e*x**5/5 + 3*a*d*e**2*x**7/7 + a*e**3* 
x**9/9 + b*d**3*x**3*asin(c*x)/3 + 3*b*d**2*e*x**5*asin(c*x)/5 + 3*b*d*e** 
2*x**7*asin(c*x)/7 + b*e**3*x**9*asin(c*x)/9 + b*d**3*x**2*sqrt(-c**2*x**2 
 + 1)/(9*c) + 3*b*d**2*e*x**4*sqrt(-c**2*x**2 + 1)/(25*c) + 3*b*d*e**2*x** 
6*sqrt(-c**2*x**2 + 1)/(49*c) + b*e**3*x**8*sqrt(-c**2*x**2 + 1)/(81*c) + 
2*b*d**3*sqrt(-c**2*x**2 + 1)/(9*c**3) + 4*b*d**2*e*x**2*sqrt(-c**2*x**2 + 
 1)/(25*c**3) + 18*b*d*e**2*x**4*sqrt(-c**2*x**2 + 1)/(245*c**3) + 8*b*e** 
3*x**6*sqrt(-c**2*x**2 + 1)/(567*c**3) + 8*b*d**2*e*sqrt(-c**2*x**2 + 1)/( 
25*c**5) + 24*b*d*e**2*x**2*sqrt(-c**2*x**2 + 1)/(245*c**5) + 16*b*e**3*x* 
*4*sqrt(-c**2*x**2 + 1)/(945*c**5) + 48*b*d*e**2*sqrt(-c**2*x**2 + 1)/(245 
*c**7) + 64*b*e**3*x**2*sqrt(-c**2*x**2 + 1)/(2835*c**7) + 128*b*e**3*sqrt 
(-c**2*x**2 + 1)/(2835*c**9), Ne(c, 0)), (a*(d**3*x**3/3 + 3*d**2*e*x**5/5 
 + 3*d*e**2*x**7/7 + e**3*x**9/9), True))
 
3.7.16.7 Maxima [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 384, normalized size of antiderivative = 1.34 \[ \int x^2 \left (d+e x^2\right )^3 (a+b \arcsin (c x)) \, dx=\frac {1}{9} \, a e^{3} x^{9} + \frac {3}{7} \, a d e^{2} x^{7} + \frac {3}{5} \, a d^{2} e x^{5} + \frac {1}{3} \, a d^{3} x^{3} + \frac {1}{9} \, {\left (3 \, x^{3} \arcsin \left (c x\right ) + c {\left (\frac {\sqrt {-c^{2} x^{2} + 1} x^{2}}{c^{2}} + \frac {2 \, \sqrt {-c^{2} x^{2} + 1}}{c^{4}}\right )}\right )} b d^{3} + \frac {1}{25} \, {\left (15 \, x^{5} \arcsin \left (c x\right ) + {\left (\frac {3 \, \sqrt {-c^{2} x^{2} + 1} x^{4}}{c^{2}} + \frac {4 \, \sqrt {-c^{2} x^{2} + 1} x^{2}}{c^{4}} + \frac {8 \, \sqrt {-c^{2} x^{2} + 1}}{c^{6}}\right )} c\right )} b d^{2} e + \frac {3}{245} \, {\left (35 \, x^{7} \arcsin \left (c x\right ) + {\left (\frac {5 \, \sqrt {-c^{2} x^{2} + 1} x^{6}}{c^{2}} + \frac {6 \, \sqrt {-c^{2} x^{2} + 1} x^{4}}{c^{4}} + \frac {8 \, \sqrt {-c^{2} x^{2} + 1} x^{2}}{c^{6}} + \frac {16 \, \sqrt {-c^{2} x^{2} + 1}}{c^{8}}\right )} c\right )} b d e^{2} + \frac {1}{2835} \, {\left (315 \, x^{9} \arcsin \left (c x\right ) + {\left (\frac {35 \, \sqrt {-c^{2} x^{2} + 1} x^{8}}{c^{2}} + \frac {40 \, \sqrt {-c^{2} x^{2} + 1} x^{6}}{c^{4}} + \frac {48 \, \sqrt {-c^{2} x^{2} + 1} x^{4}}{c^{6}} + \frac {64 \, \sqrt {-c^{2} x^{2} + 1} x^{2}}{c^{8}} + \frac {128 \, \sqrt {-c^{2} x^{2} + 1}}{c^{10}}\right )} c\right )} b e^{3} \]

input
integrate(x^2*(e*x^2+d)^3*(a+b*arcsin(c*x)),x, algorithm="maxima")
 
output
1/9*a*e^3*x^9 + 3/7*a*d*e^2*x^7 + 3/5*a*d^2*e*x^5 + 1/3*a*d^3*x^3 + 1/9*(3 
*x^3*arcsin(c*x) + c*(sqrt(-c^2*x^2 + 1)*x^2/c^2 + 2*sqrt(-c^2*x^2 + 1)/c^ 
4))*b*d^3 + 1/25*(15*x^5*arcsin(c*x) + (3*sqrt(-c^2*x^2 + 1)*x^4/c^2 + 4*s 
qrt(-c^2*x^2 + 1)*x^2/c^4 + 8*sqrt(-c^2*x^2 + 1)/c^6)*c)*b*d^2*e + 3/245*( 
35*x^7*arcsin(c*x) + (5*sqrt(-c^2*x^2 + 1)*x^6/c^2 + 6*sqrt(-c^2*x^2 + 1)* 
x^4/c^4 + 8*sqrt(-c^2*x^2 + 1)*x^2/c^6 + 16*sqrt(-c^2*x^2 + 1)/c^8)*c)*b*d 
*e^2 + 1/2835*(315*x^9*arcsin(c*x) + (35*sqrt(-c^2*x^2 + 1)*x^8/c^2 + 40*s 
qrt(-c^2*x^2 + 1)*x^6/c^4 + 48*sqrt(-c^2*x^2 + 1)*x^4/c^6 + 64*sqrt(-c^2*x 
^2 + 1)*x^2/c^8 + 128*sqrt(-c^2*x^2 + 1)/c^10)*c)*b*e^3
 
3.7.16.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 711 vs. \(2 (259) = 518\).

Time = 0.30 (sec) , antiderivative size = 711, normalized size of antiderivative = 2.48 \[ \int x^2 \left (d+e x^2\right )^3 (a+b \arcsin (c x)) \, dx=\frac {1}{9} \, a e^{3} x^{9} + \frac {3}{7} \, a d e^{2} x^{7} + \frac {3}{5} \, a d^{2} e x^{5} + \frac {1}{3} \, a d^{3} x^{3} + \frac {{\left (c^{2} x^{2} - 1\right )} b d^{3} x \arcsin \left (c x\right )}{3 \, c^{2}} + \frac {b d^{3} x \arcsin \left (c x\right )}{3 \, c^{2}} + \frac {3 \, {\left (c^{2} x^{2} - 1\right )}^{2} b d^{2} e x \arcsin \left (c x\right )}{5 \, c^{4}} + \frac {6 \, {\left (c^{2} x^{2} - 1\right )} b d^{2} e x \arcsin \left (c x\right )}{5 \, c^{4}} + \frac {3 \, {\left (c^{2} x^{2} - 1\right )}^{3} b d e^{2} x \arcsin \left (c x\right )}{7 \, c^{6}} - \frac {{\left (-c^{2} x^{2} + 1\right )}^{\frac {3}{2}} b d^{3}}{9 \, c^{3}} + \frac {3 \, b d^{2} e x \arcsin \left (c x\right )}{5 \, c^{4}} + \frac {9 \, {\left (c^{2} x^{2} - 1\right )}^{2} b d e^{2} x \arcsin \left (c x\right )}{7 \, c^{6}} + \frac {{\left (c^{2} x^{2} - 1\right )}^{4} b e^{3} x \arcsin \left (c x\right )}{9 \, c^{8}} + \frac {\sqrt {-c^{2} x^{2} + 1} b d^{3}}{3 \, c^{3}} + \frac {3 \, {\left (c^{2} x^{2} - 1\right )}^{2} \sqrt {-c^{2} x^{2} + 1} b d^{2} e}{25 \, c^{5}} + \frac {9 \, {\left (c^{2} x^{2} - 1\right )} b d e^{2} x \arcsin \left (c x\right )}{7 \, c^{6}} + \frac {4 \, {\left (c^{2} x^{2} - 1\right )}^{3} b e^{3} x \arcsin \left (c x\right )}{9 \, c^{8}} - \frac {2 \, {\left (-c^{2} x^{2} + 1\right )}^{\frac {3}{2}} b d^{2} e}{5 \, c^{5}} + \frac {3 \, {\left (c^{2} x^{2} - 1\right )}^{3} \sqrt {-c^{2} x^{2} + 1} b d e^{2}}{49 \, c^{7}} + \frac {3 \, b d e^{2} x \arcsin \left (c x\right )}{7 \, c^{6}} + \frac {2 \, {\left (c^{2} x^{2} - 1\right )}^{2} b e^{3} x \arcsin \left (c x\right )}{3 \, c^{8}} + \frac {3 \, \sqrt {-c^{2} x^{2} + 1} b d^{2} e}{5 \, c^{5}} + \frac {9 \, {\left (c^{2} x^{2} - 1\right )}^{2} \sqrt {-c^{2} x^{2} + 1} b d e^{2}}{35 \, c^{7}} + \frac {{\left (c^{2} x^{2} - 1\right )}^{4} \sqrt {-c^{2} x^{2} + 1} b e^{3}}{81 \, c^{9}} + \frac {4 \, {\left (c^{2} x^{2} - 1\right )} b e^{3} x \arcsin \left (c x\right )}{9 \, c^{8}} - \frac {3 \, {\left (-c^{2} x^{2} + 1\right )}^{\frac {3}{2}} b d e^{2}}{7 \, c^{7}} + \frac {4 \, {\left (c^{2} x^{2} - 1\right )}^{3} \sqrt {-c^{2} x^{2} + 1} b e^{3}}{63 \, c^{9}} + \frac {b e^{3} x \arcsin \left (c x\right )}{9 \, c^{8}} + \frac {3 \, \sqrt {-c^{2} x^{2} + 1} b d e^{2}}{7 \, c^{7}} + \frac {2 \, {\left (c^{2} x^{2} - 1\right )}^{2} \sqrt {-c^{2} x^{2} + 1} b e^{3}}{15 \, c^{9}} - \frac {4 \, {\left (-c^{2} x^{2} + 1\right )}^{\frac {3}{2}} b e^{3}}{27 \, c^{9}} + \frac {\sqrt {-c^{2} x^{2} + 1} b e^{3}}{9 \, c^{9}} \]

input
integrate(x^2*(e*x^2+d)^3*(a+b*arcsin(c*x)),x, algorithm="giac")
 
output
1/9*a*e^3*x^9 + 3/7*a*d*e^2*x^7 + 3/5*a*d^2*e*x^5 + 1/3*a*d^3*x^3 + 1/3*(c 
^2*x^2 - 1)*b*d^3*x*arcsin(c*x)/c^2 + 1/3*b*d^3*x*arcsin(c*x)/c^2 + 3/5*(c 
^2*x^2 - 1)^2*b*d^2*e*x*arcsin(c*x)/c^4 + 6/5*(c^2*x^2 - 1)*b*d^2*e*x*arcs 
in(c*x)/c^4 + 3/7*(c^2*x^2 - 1)^3*b*d*e^2*x*arcsin(c*x)/c^6 - 1/9*(-c^2*x^ 
2 + 1)^(3/2)*b*d^3/c^3 + 3/5*b*d^2*e*x*arcsin(c*x)/c^4 + 9/7*(c^2*x^2 - 1) 
^2*b*d*e^2*x*arcsin(c*x)/c^6 + 1/9*(c^2*x^2 - 1)^4*b*e^3*x*arcsin(c*x)/c^8 
 + 1/3*sqrt(-c^2*x^2 + 1)*b*d^3/c^3 + 3/25*(c^2*x^2 - 1)^2*sqrt(-c^2*x^2 + 
 1)*b*d^2*e/c^5 + 9/7*(c^2*x^2 - 1)*b*d*e^2*x*arcsin(c*x)/c^6 + 4/9*(c^2*x 
^2 - 1)^3*b*e^3*x*arcsin(c*x)/c^8 - 2/5*(-c^2*x^2 + 1)^(3/2)*b*d^2*e/c^5 + 
 3/49*(c^2*x^2 - 1)^3*sqrt(-c^2*x^2 + 1)*b*d*e^2/c^7 + 3/7*b*d*e^2*x*arcsi 
n(c*x)/c^6 + 2/3*(c^2*x^2 - 1)^2*b*e^3*x*arcsin(c*x)/c^8 + 3/5*sqrt(-c^2*x 
^2 + 1)*b*d^2*e/c^5 + 9/35*(c^2*x^2 - 1)^2*sqrt(-c^2*x^2 + 1)*b*d*e^2/c^7 
+ 1/81*(c^2*x^2 - 1)^4*sqrt(-c^2*x^2 + 1)*b*e^3/c^9 + 4/9*(c^2*x^2 - 1)*b* 
e^3*x*arcsin(c*x)/c^8 - 3/7*(-c^2*x^2 + 1)^(3/2)*b*d*e^2/c^7 + 4/63*(c^2*x 
^2 - 1)^3*sqrt(-c^2*x^2 + 1)*b*e^3/c^9 + 1/9*b*e^3*x*arcsin(c*x)/c^8 + 3/7 
*sqrt(-c^2*x^2 + 1)*b*d*e^2/c^7 + 2/15*(c^2*x^2 - 1)^2*sqrt(-c^2*x^2 + 1)* 
b*e^3/c^9 - 4/27*(-c^2*x^2 + 1)^(3/2)*b*e^3/c^9 + 1/9*sqrt(-c^2*x^2 + 1)*b 
*e^3/c^9
 
3.7.16.9 Mupad [F(-1)]

Timed out. \[ \int x^2 \left (d+e x^2\right )^3 (a+b \arcsin (c x)) \, dx=\int x^2\,\left (a+b\,\mathrm {asin}\left (c\,x\right )\right )\,{\left (e\,x^2+d\right )}^3 \,d x \]

input
int(x^2*(a + b*asin(c*x))*(d + e*x^2)^3,x)
 
output
int(x^2*(a + b*asin(c*x))*(d + e*x^2)^3, x)